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We argue for a recently observed puzzling multiferroic behavior in s=l/2 ID chain cuprate 
LiCu2 02 with edge-shared arrangement of Cu04 plaquettes and incommensurate spiral spin or- 
dering can be consistently explained if one takes into account the nonrelativistic exchange-induced 
electric polarization on the Cu^"*" centers substituting for the positions native for the Cu^"'"-ions. 
These substituent centers are proved to be an effective probe of the spin incommensurability and 
magnetic field effects. 



Recent observations of multiferroic behaviour concomi- 
tant the incommensurate spin spiral ordering in chain 
cuprates LiCuV04^'^'^ and LiCu202^ challenge the mul- 
tiferroic community. At first sight, these cuprates seem 
to be prototypical examples of ID spiral-magnetic ferro- 
electrics revealing the relativistic mechanism of "ferro- 
electricity caused by spin-currents"^ with the textbook 
expression for the uniform polarization induced by a 
spin spiral with the wave vector Q: P oc [es x Q] , where 
63 is a vector orthogonal to the spin spiral plane^ or 
Pij oc [Rij X [Si X Sj]], where denotes the vector 
connecting the two sites and [Si x Sj] is a local spin 
currenti^ However, the both systems reveal a mysteri- 
ous behavior with conflicting results obtained by differ- 
ent groups. Indeed, Yasui et al^ claim the LiCuV04 re- 
veals clear deviations from the predictions of spin-current 
models^'^ while Schrettle et al^ assure of its applicabil- 
ity. In contrast to LiCuV04, the LiCu202 shows up a 
behavior which is obviously counterintuitive within the 
framework of spiral-magnetic ferroelectricityi^ It is worth 
noting that at variance with Park et alA, Naito et al^ 
have not found any evidence for ferroelectric anoma- 
lies in LiCu202. Such a discrepancy one observes in 
microscopic model approaches as well. The relativis- 
tic LSDA calculations^ seemingly explain the LiCuV04 
data'^ but fail in case of LiCu202. However, a detailed 
analysis of relativistic effects for the system of gg-holes 
in a perfect chain structure of edge-shared CUO4 pla- 
quettes as in LiCuV04 and LiCu202 shows that the in- 
chain spin current does not produce an electric polar- 
ization because of an exact cancellation of two Cu-0- 
Cu paths4 Moreover, recently we have shown^ that the 
multiferroicity in LiCuV04 may have nothing to do with 
relativistic effects and can be consistently explained, if 
the nonrelativistic exchange-induced electric polariza- 
tion on the out-of-chain Cu^^ centers substituting for 
Li- ions in LiVCu04 is taken into account. Below we ar- 
gue that a similar mechanism which takes into account 
the exchange-induced electric polarization on the Cu^+ 
centers, substituting unexpectedly for Cu^+-ions, is at 
work in LiCu202. 

LiCu202 is orthorhombic mixed- valent compound with 
copper ions in the Cu^+ and Cu^+ valence statesplS The 



unit cell contains four magnetic Cu ions belonging to 
two pairs of Cu02 chains formed by edge-shared Cu^+04 
plaquettes running along the crystallographic 6-axis and 
linked by the LiO^ double chains. Alternating double 
parallel chains, containing either Li or Cu atoms, form 
the sheets which are interconnected by Cu^+ in 0-Cu-O 
dumbbells. 

The first experimental evidence of magnetic incom- 
mensurability in LiCu202 was obtained independently 
by Gippius et al^ and Masuda et alP from •^■'^Li NMR 
and neutron diffraction measurements, respectively. Any 
spins related by a translation along the c axis and a 
axis are parallel and antiparallel to each other, respec- 
tively. A good fit to neutron diffraction data was ob- 
tained with all spins confined to the ab crystallographic 
plane^'^ thus forming a6-plane spin spirals running along 
b axis:S(i/) = S'(cos6', sin0, 0), where = qy + a, a is 
a phase shift. Park et al."^ have found that the incom- 
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FIG. I: Direction of ferroelectric polarization in LiCu2 02 for 
different spin spiral plane orientation. 



mensurate spin ordering in LiCu202 below T^r « 23 K 
is accompanied by a ferroelectric transition with a puz- 
zling anisotropy and field dependence which are repro- 
duced schematically in Fig. 1. First of all, the electric 
polarization in zero field is directed along the c axis im- 
plying in accordance with the concept of spin current 
induced ferroelectricity that the spiral spins lie in the 
6c-planc in sharp contrast with earlier neutron diffrac- 
tion data.^^ When a magnetic field applied along the b 
axis (see Refi for the making use of a, 6 notations in 
ab-twinned crystal), Pc decreases and Pa increases, im- 
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FIG. 2: An idealized view of crystal structure of 
LiCu202(upper bilayer). A "left" site impurity center with 
Cu^^-ion substituted for Cu^^-ion is inbetween upper and 
lower Cu02 chains from the same unit cell. Shown is the hole 
density distribution in dyz orbital. The exchange induced 
dipole moments are shown by arrows. 



plying that the Cu^+ spin spiral plane flips from the be 
to ab plane, resulting in a flip of the polarization from 
the c to a axis. It is expected that h || c may flip the 
spiral plane from the be to ab plane, so that P may flip 
from the c to a axis with he- However, this is com- 
pletely in contrast with the observations'* that he en- 
hances Pc and hj, is the one inducing the P flip from 
the c to a axis. The appearance of Pa with ha is also 
counterintuitive within the framework of the relativis- 
tic spiral-magnetic ferroelectricityi^i^ These unexpected 
magnetic field effects raise doubts about the validity of 
the scenario of relativistic spin curent spiral-magnetic fer- 
roelectricity and point to another, probably the out-of- 
spin-chain origin of the magnetoelectric coupling. In this 
connection it is worth noting that the thermogravimet- 
ric analysis revealed the LiCu202 samples had a lower 
content of Cu ions than follows from the stoichiomet- 
ric formulaii^ Chemical disorder and a Cu deficiency by 
asmuch as x = 16% are inherently present. The "sur- 
plus" Li+ ions in LiCu202 occupy Cu^+ sites, due to a 
good match of ionic radii (0.68 and 0.69 A, respectively). 
The charge compensation requires that the introduction 
of nonmagnetic Li+ ions into the double chains is ac- 
companied by a transfer of the S=l/2-carrying Cu^+ 
ions onto the Cu"*" interchain sites »^ At first sight it 
seems improbable because of different coordination pref- 
erences. However, the actual coordination of the native 
Cu"*" interchain site approaches most likely to an axially 
distorted square, or rhombic coordination due to an ex- 
tremely small inter-dumbbell separation {d « 2. 86 A) as 
compared with other 0-Cu'^-O dumbbell bearing com- 
pounds (e.g. YBa2Cu06, d ~ 3.8A)^ In other words, 
the Zhang-Rice singlet within the Cu02 chains becomes 
unstable with respect to a hole transfer to one of neigh- 
boring Cu''" sites. Details of this instability will be dis- 
cussed elsewhere. If the doped hole would be remain 



in the Cu02 chains, dimer-type effects as in other hole 
doped chains would be observed experimentally. Also the 
spiral state observed in the neutron diffraction would be 
strongly disturbed by the presence of these holes. What 
is the ground state of the single hole configuration of 
Cu^+ ion in the native Cu"*" interchain sites? Purely elec- 
trostatic arguments made within the framework of the 
point charge model, supported by account for Cu 3d-0 2p 
covalency, point to a competition of and dyz orbitals 
while strong intra-atomic s-dz2 hybridization singles out 
the dyz orbital to be a main candidate for the ground 
state. The Cu^+ substituents in native Cu"'^"'" positions 
form strongly polarizable entities which electric polar- 
ization due to a parity-breaking exchange interactioni^ 
with Cu^+ spin spirals explains all the puzzles observed 
by Park et al."^ This unconventional exchange coupling 
can be easily illustrated for, e.g., the one-particle (elec- 
tron/hole) center in a crystallographically centrosymmet- 
ric position of a magnetic crystal when all the particle 
states can be of definite spatial parity, even (g) or odd 
(u), respectively. Having in mind the 3d centers we'll 
assume an even-parity ground state \g). For simplicity 
we restrict ourselves by only one excited odd-parity state 
\u). The exchange coupling with the surrounding spins 
can be written as follows: 

v;. = 5]/(R„)(s.s„), (1) 

n 

where /(Rn) is an orbital operator with a matrix 

/(R„)- f^^«{5-). l'"!SiV (2) 

The parity-breaking off-diagonal part of the exchange 
coupling can lift the center of symmetry and mix \g) and 
\u) states giving rise to a nonzero electric dipole polar- 
ization of the ground state 

P = 2cg„(g|er|u) =^n„(s-S„) (3) 

n 

with n„ = 2/g„(R„)(g|er|u)/A„g (A„g = e„ - e^). 

Strictly speaking, the parity-breaking exchange cou- 
pling of native Cu^"*" center in Cu02 chain (hole ground 
state \g) cx dj;y) with neighbouring Cu^^ substituent 
(hole ground state \g) cx dyz) will result in the a6-plane 
electric polarization of CUO4 chain plaquettes and the 
c-axis polarization of the Cu^+ substituent. 

Unit cell of LiCu202 contains two types (left and right) 
of native Cu^+-positions (see Fig.[T]) with four neigh- 
bouring Cu^+ centers in the two CUO2 chains. Within 
the framework of our model the both "left" A- type and 
"right" B-type substituent centers differ by the spin spi- 
ral phase shift a = 7r/2 and a — — 7r/2 with respect 
to the lower chain, and by orientation of the generated 
electric dipole moments: da{A) = —da{B) = (i/\/2, 
4(A) = di,{B) = d/V2, dc{A) = 4(B) = 4. Here we ig- 
nore the weak infiuence of the adjacent chain in the third 
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Cu02 chain belonging to the adjacent bilayer. According 
to LDA calculation there is practically no hybridization 
with that chain.— 

To describe different configurations of the spin neigh- 
bourhood for a Cu^+ substituent (see Fig. [2]) we intro- 
duce four basic vectors similarly to conventional ferro- 
and antiferromagnetic vectors as follows: 

F{y) = [Si + S2 + S3 + S4]; G{y) = [Si - S2 -f S3 - S4]; 



A{y) = [Si + S2 - S3 - S4]; C{y) = [Si - S2 - S3 + S4] 

with a kinematic constraint: (F • A) = (C • G)=0 valid 
for two identical spirals irrespective of their phase shift. 
Then the electric polarization induced by the parity- 
breaking exchange coupling of Cu^"*" substituent with a 
complete set of four neighboring in-chain Cu^+ ions 1-4 
(see Fig.[T]) can be written as follows: 



Pa 



(s- A); Pfc -4(s-C) 



Spin polarization of Cu^+ substituent spin can be eas- 
ily found within the framework of a weak coupling ap- 
proximation, if one take the most general form of the 
impurity-spiral ground state {gg) exchange interaction 



VsS = 



J2 si (i) S, = (s • Ho) 



(4) 



i=l-4 



where Hq is an effective magnetic field, acting on the 

Cu^+ substituent, Iaa{i) = Iaa,Ixz{i) = Ixz;Ixy{l) = 
^xyi'^) ~ -^xyi^) — -^xyi^) — -^xy'i zy{^) ~ ^zy^^ 

I .' i = —Izyi^) = Izy are a symmetric matrix of the 
exchange integrals. Thus for the effective field we obtain 



Ho(2/) = IfF+ IgG 



(5) 
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We start with the a6-plane spiral ordering of Cu^+ 
spins in the Cu02 chains of LiCu202 , which deduced 
from neutron diffraction data in zero external field^- and 
assume T = 0. For zero external magnetic field or for 
a field directed along the c axis and for a — ±7r/2 the 
electric polarization of the y-th Cu^"*" substituent center 
oscillates as follows 



Pciy) = 



H{y) 



-[{Ixx - Iyy)ucos{2qy) ± 2IxyV ain{2qy)] , 

(6) 

where u — cos{^),v — sin(^), and for h = 

H{y) = 2V25[(4\ + 4 + Il,)u' + (2^ + ll^) 



Tmxx+Iyy)iIxx-Iyy)+llW-I%v')smi2qy)]-^ . (7) 

First, it should be noted that the both "left" A- type 
{a — 7r/2) and "right" B-type {a = —tt/2) substituent 
positions contribute equally to a macroscopic polariza- 
tion Pc- On the other hand, it means that Pa vanishes 
due to an exact compensation of A-type and B-type con- 
tributions since da{A) — —da{B). For we arrive at a 
strict cancellation of the net electric polarization given 
a ±7r/2 due to opposite signs of the antisymmetric 
part of the effective field. Moreover, this cancellation 
hold itself also under an external magnetic field irrespec- 
tive of its direction. Second, we note that a nonzero 
electric polarization for the substituent center 12-Cu^^- 
34 can be related only with the anisotropic substituent- 
spiral exchange coupling. The net polarization (Pdy)) 
seems to be rather weak because of several reduction ef- 
fects: i) the existence of non-compensated non-oscillatory 
contribution of isotropic exchange to the effective mag- 
netic field ([7]); ii) a quadratic or cubic dependence of 
{Pc{y)) on the exchange anisotropy parameters. In or- 
der to demonstrate the role of the anisotropic exchange 
we adopt a relation between the anisotropy parameters 
predicted by a simple nearest-neighbor magneto-dipole 
model: (hx-Iyy) = 0, (I^x = hy) = V2Ixy The depen- 
dence of {Pc{y)) on the ratio 5 = Ixyjlxx appears to be 
strongly nonlinear, being approximately cx for a small 
anisotropy. Only a strong anisotropy 8 ^ \ provides the 
magnitudes of {Pc{y)) comparable with that of {Pa{y)) 
in LiCuV04. The typical field dependence of {Pc{y)) is 
shown in Fig. [3] given qb/2 = 0.1727r which corresponds 
to a pitch angle « 62°. i^^^ 

A magnetic field h || a induces in LiCu202 a ab-bc 
spin-flop transition to the phase with a 6c-plane spiral 
ordering. Interestingly, that irrespective of the fleld di- 
rection a 6c-plane spin spiral ordering, similarly to that 
of a6-plane one, supports only a c-axis orientation of both 
local and net electric polarizations, which expressions 
can be easily obtained from their a6-axis counterparts. 
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FIG. 3: The field dependence of {Pc{y)) (in units of dc) for 
a&-plane spiral: Ixy/ Ixx = 0.4 
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<P>: ac-plane spiral 




FIG. 4: The field dependence of {Pa.bijj)) (in units of d) for 
h II b {Ixy/Ixx = 0.4) for the ac-plane spiral. 



if one makes the interchange — > /iz, Ixy Izy 

It 

is worth noting that at variance with the a6-plane spin 
spiral ordering the c-axis orientation of net electric po- 
larization for 6c-plane spin arrangement agrees with the 
predictions of the spin current scenario. Thus both the 
local P(i/) and the averaged electric polarizations (P(i/)) 
for ab and be plane spin spirals lie along the c-axis even 
in zero magnetic field. It's quite another matter for the 
ac-plane spin spiral arrangement which can be a result 
of a spin-flop transition in an external magnetic field 
directed along 6-axis. Contributions of the A- and B- 
type centers to Pa and Pc are strictly opposite in sign, 
that means their cancellation for Pc and doubling for Pa . 
From the other hand, for the first time, the Pf, compo- 
nent of electric polarization appears to be nonzero. Thus, 
in contrast with two preceding instances both the local 
P(j/) and the averaged electric polarizations (P(y)) for 
ac-plane spin spirals lie in ab plane even in zero mag- 
netic field. Moreover, we arrive at a simple relation be- 
tween the a- and the 6-components of the electric po- 
larization: Pb/Pa ~ —v/u = —tan{qb/2), that corre- 
sponds to Pb/Pa ~ —0.6 given the pitch angle qb « 62°. 
Fig. m shows the field dependence of {Pa,biy)) for h || b 
{Ixy/Ixx = 0.4) for the ac-plane spiral. 

The mechanism of impure ferroelectricity we discuss 
does consistently explain all the puzzles of the magneto- 
electric effect observed in LiCu202 by Park et aZ.^(see 



Fig. 1). First of all the model explains the c-axis direc- 
tion of the spontaneous electric polarization emerging be- 
low the spiral-magnetic ordering temperature within the 
framework of a dominant ab plane Cu^^ spin arrange- 
ment, proposed earlier from neutron diffraction dataJ^ 
We argue that an external field h || b induces spin- flop 
transition with the Cu^+ spin spiral plane flipping from 
the ab to the ac plane accompanied by the flipping of net 
electric polarization P from the c-axis to the a6-plane 
where the relation imbetween b- and a-components is de- 
termined by the actual pitch angle. The twin structure 
observed in ab plane of the LiCu202 crystal and fer- 
roelectric domain effects^ make the field dependence of 
electric polarization quite complex. Indeed, an external 
field h II b induces different spin- flop transitions in differ- 
ent twins: ab—^ac and ab-^bc, respectively. Only in the 
former twins we deal with Pc-^Pab flipping of ferroelec- 
tric moment, while in the latter twins the polarization 
remains oriented along the c-axis, though having a var- 
ied magnitude as compared with the ab plane spin spiral. 
Such a behavior is observed in experiments by Park et 
al.,^ with a relation between the in-plane components of 
polarization which is close to a theoretically predicted 
value 0.6. At variance with LiCuV04 the spontaneous 
{h = 0) electric polarization in the low-temperature 
spiral phase of LiCu2 02 does depend not only on the 
pitch angle (qb) and the relation inbetween the values 
of exchange anisotropy parameters, but also on the rela- 
tive magnitude of exchange anisotropy as compared with 
isotropic exchange. Namely this feature is believed to de- 
termine the relatively small magnitude of the multiferroic 
effect in LiCu202 as compared with LiCuV04ji 

Thus we conclude that at variance with the relativistic 
spin current model the nonrelativistic parity breaking 
exchange induced polarization for the centers formed 
by Cu^"*" substituted for Cu^+ in nonstoichiometric 
LiCu202with a simple zero-field afe-plane spiral order- 
ing can be a natural electronic source of multiferroicity 
found by Park et al.^ in this cuprate. 
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